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ABSTRACT 

We present evidence which indicates the luminosity functions of star forming regions, 
open star clusters and the field dips between spectral types M7-M8. We attribute this 
to a sharp local drop in the luminosity-mass relation and speculate that this is caused 
by the beginning of dust formation in the atmospheres of objects in this effective 
temperature regime. This effect is not predicted by the current generation of low mass 
stellar/substellar evolutionary models. If our interpretation is correct then this result 
has important implications for investigations concerned with the mass functions of 
star forming regions and young open star clusters. For example, it suggests that some 
brown dwarfs have higher masses than previously thought. 

Key words: stars: late-type, low-mass, brown dwarfs, atmospheres - infrared: stars. 



1 INTRODUCTION 

Over the last few years deep far-red optical and infrared sur- 
veys of young open clusters and star forming regions have 
revealed large populations of brown dwarfs (e.g Bouvier et 
al. 1998, Pinfield et al. 2000, Magazzu et al. 1998, Luhman 
et al. 1999, 2000, Lodieu et al. 2002). Indeed, several of 
the most recent investigations may have unearthed signif- 
icant numbers of free-floating objects with masses less than 
the deuterium burning limit (M< O.OI2M0; e.g ct— Orionis, 
Zapatero-Osorio et al. 2000, ^— Orionis, Lucas & Roche 
2000, Muench et al. 2002), some possibly having masses as 
low as M~ O.OO6M0 (Lucas & Roche 2001). The form of the 
initial mass function (IMF; dN/dm) in star forming regions, 
young open clusters and in the field, across and below the 
stellar/substellar boundary and down to at least O.O3M0, is 
found to be well matched by a slowly rising powerlaw model 
with an index in the range a ~ — 1.5 (e.g Bouvier et al. 

1998, Luhman et al. 1999, Hambly et al. 1999, Reid et al. 

1999, Bejar et al. 2001, Bobbie et al. 2002). 

However, there is now some evidence to suggest that 
below M~ O.O3M0 the IMF departs from this simple de- 
scription. Two independent studies of the 0— Orionis star 
forming region appear to show a significant secondary peak 
in the IMF close to the deuterium burning limit (Muench 
et al. 2002, Lucus & Roche 2000). The IMF of a-Orionis 
may also rise steeply at these masses (Bejar et al. 2001). We 
note that standard star formation theories (e.g Boss 1998, 
Low & Lynden-Bell 1976) predict that the minimum mass 
of gas which is unstable against collapse to stellar densities 
is of the order O.OlMg (e.g Silk 1977). Muench et al. (2002) 
have suggested that another formation mechanism may con- 



tribute to the IMF at these very low masses to produce the 
observed secondary peak. 

Of course, derivations of the IMF are dependent on the 
robustness of the conversion from the observational to the 
theoretical plane. Studies focusing on the very youngest star 
clusters and star forming regions are undoubtedly compli- 
cated by residual dust and gas along the lines of sight and 
the presence of T-Tau disks (e.g Lucas & Roche 2000). There 
are also uncertainties in the theoretical models, particularly 
at these youngest ages. For example, at ages of less than a 
few Myrs, predicted TeffS and luminosities are found to be 
highly sensitive to the choice of initial conditions and the 
adopted treatment of convection (e.g Baraffe et al. 2002). 
Furthermore, there is now little doubt that the spectral en- 
ergy distributions of the coolest stellar and substellar ob- 
jects are strongly influenced by the presence of dust grains 
in their atmospheres (e.g Tsuji et al. 1996, Jones & Tsuji 
1997, Chabrier et al. 2000, Marley et al. 2002). While the 
observed trends in colours of the late-M dwarfs through the 
L dwarfs to the T dwarfs can now be qualitatively repro- 
duced by the latest model atmospheres which incorporate a 
treatment for both the condensation and sedimentation of 
species such as AI2O3, MgSiOa and Fe (e.g Tsuji et al. 2001), 
the theory of dust grain formation and gravitational settling 
in stellar photospheres is still very much in its infancy. The 
effects of these processes on the evolution of substellar ob- 
jects have not yet been investigated in a fully self-consistent 
manner. 

In the next section of this work we present evidence 
for a dip in the luminosity function (LF) between spectral 
types M7-M8. In Section 3 we attribute this to a sharp lo- 
cal change in the form of the luminosity-mass relation. We 
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Figure 1. A subsection o{ 1,1 — K plot for the Pleiades pre- 
sented in Jameson et al. (2002). Candidate low mass members 
with photometry consistent with being single and binary objects 
are shown by the filled and open circles respectively. The 125Myr 
NEXTGEN isochrone with discrete mass points marked (crosses) 
and a corresponding sequence for equal mass binaries are also 
shown (dashed and dotted lines respectively). The sparsely pop- 
ulated region of the cluster sequence between 3.5 < I — K < 3.7 
is highlighted (thick grey line). 

construct semi-empirical magnitude-mass relations for the 
Pleiades and a— Orionis clusters using the available obser- 
vational data. In Section 4 we speculate that the change is 
due to the onset of dust formation in the atmospheres of 
cool stellar and substellar objects and discuss the implica- 
tions of our result for current investigations concerned with 
determining the form of the IMF of star forming regions and 
young open clusters. It should be noted that throughout this 
work, we refer to 7c as /. 

2 EVIDENCE FOR THE MISSING M DWARFS 
2.1 The Pleiades 

In our recent work concerned with placing stringent con- 
straints on the core radius of the substellar members of the 
Pleiades and on the shape of the cluster mass function across 
and below the stellar/substellar boundary we have compiled 
a large, complete, magnitude limited sample of candidate 
brown dwarfs (Jameson et al. 2002). As part of this effort 
we have constructed an 1,1 — K colour-magnitude diagram 
(CMD) using 48 candidate substellar Pleiads which has re- 
vealed an apparent sparsity of cluster members with colours 
3.5 < I — K < 3.7 (or 2.5 ^ / — J ^ 2.7 from a comparison 
with the colours of field dwarfs). Using data from Pinfield 
et al. (2002, in prep) and Tables 2 and 3 of Jameson et 
al. (2002) we replot the relevant section of the CMD in Fig- 
ure 1 where candidate single and unresolved binary low mass 
members are shown as filled and open circles respectively. 
The location of the 125Myr NEXTGEN isochrone (Baraffe 
et al. 1998) for solar metalicity and the corresponding equal 
mass binary track, after correcting for extinction, redden- 
ing and a distance of A/ — 0.07, E(7 — K)= 0.06 and (m- 
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Figure 2. The /-band LF, in the interval 17.7 < / < 19.5, of 
the candidate Pleiads (excluding likely binaries) from the 7.6 sq. 
degree region of the cluster included in the Jameson et al. (2002) 
study. A powerlaw model of the LF is also shown (heavy dashed 
line; see text for details). We find the discrepancy between the 
observed and model LF to be significant at a level of confidence 
^ 90%. 

M)o = 5.53 (Crawford & Perry 1976) respectively are also 
shown (dashed and dotted lines respectively; see arguments 
in Bobbie et al. 2002 for our choice of cluster distance). It 
can be seen that the region of the cluster sequence which we 
find to be sparsely populated (highlighted by the gray shad- 
ing) spans ~ 0.5 magnitudes at /. Available infrared data 
(Jameson et al. 2002; Pinfield et al. 2002, in prep.) reveal 
it to extend over ~ 0.3 magnitudes at J, H and K respec- 
tively, suggesting it spans ~ 0.3 magnitudes in bolometric 
luminosity. 

The /-band LF over the interval 17.7 < / < 19.50 for 
the 7.6 sq. degrees of the Pleiades included in the Jameson 
et al. (2002) study is shown in Figure 2. Note that those 
candidates lying close to the binary sequence have been ex- 
cluded from this determination (see later). Excluding the 
magnitude bin which has an obvious deficit of members, we 
have modelled log N versus / with a powerlaw (see Muench 
et al. 2002 for justification), represented in Figure 2 by the 
dashed line. A test indicates that the observed LF, in- 
cluding the bin with the deficit of members, is inconsistent 
with this model at a level of confidence ~ 90% (for the 
adopted bin width of 0.45 magnitudes). We note that the 
colour, I ~ K ^ 3.6, at the midpoint of the LF dip corre- 
sponds to a spectral type of M7-M8 (Leggett 1992) and an 
Tcff « 2700K (NEXTGEN model). 

2.2 6— and cr Orionis 

Both the cr— Orionis star cluster and the S— Orionis (Trapez- 
ium) star forming region, which lie at ~ 400pc (e.g Blaauw 
1991, Warren & Hesser 1978), are considerably younger than 
the Pleiades with estimated ages of 1.7-7 Myrs and 0.3-2 
Myrs respectively (e.g Ali & Depoy 1995, Hillenbrand 1997, 
Brown, de Geus & de Zeeuw 1994, Warren & Hesser 1978). 
Bejar et al. (2001) present the results of a 847 square arcmin 
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Figure 3. The magnitude-mass relation for the Pleiades low- 
mass stars and brown dwarfs as determined from the observations. 
The purely theoretical 125Myr NEXTGEN and 120 Myr DUSTY 
relations of Baraffe et al. (1998) and Chabrier et al. (2000) are 
also shown (solid and dashed lines respectively). The location of 
the dip in the LF is highlighted (dotted lines). 

IZJ band survey of the cr— Orionis region, which is complete 
to / = 21.5 and has unearthed 64 candidate members. Their 
Figure 1, the 1,1 — J CMD for the candidates, shows a spar- 
sity of objects with 2.6 < I ~ J < 3.0, virtually coincident 
in colour with the LF dip observed in the sequence of the 
Pleiades. Lucas & Roche (2000) present the results of an /, 
J and H band survey of Orionis complete to iif ~ 18.0, 
in which 515 sources were detected in a 33 sq. arcminute re- 
gion. The vast majority of these are believed to be members 
of the Trapezium. A J — H , I — J colour-colour plot of these 
sources, constructed using dereddened photometry, also re- 
veals a sparsity of objects in the range 2.6 < I — J < 2.8 
(see their Figure 3). Their /—band LF for the Trapezium re- 
gion is falling into the substellar regime but appears to rise 
again towards a secondary peak at M_r « 10.5 {H « 16.5, 
J « 17.5). The A'— band LF derived from an independent 
near-IR survey of the region similarly displays such a fall in 
the substellar regime, before apparently rising towards a sec- 
ondary peak at /i' « 15.5. Muench et al. (2002), who argue 
that the theoretical luminosity- mass relation is reasonably 
robust in the substellar regime, interprete this apparent sec- 
ondary peak as a rise in the IMF near the deuterium burning 
limit. 

2.3 Other clusters and the field 

Further evidence for drops in LFs in this Toff regime can 
be found elsewhere. For example, Barrado y Navascues et 
al. (2001) present the results of a VRIZ survey, complete 
to I ^ 18, of the 53Myr old open cluster IC2391. In a 2.5 
square degree area, they identify 132 probable and possible 
cluster members, although as follow-up data are available for 
only a limited number of these ~ 25% are likely to be non- 
members. From their Figure 3, the I,R — I CMD, it can be 
seen that the number of candidates drops dramatically just 
above the completeness limit of their survey, before increas- 
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Figure 4. The magnitude-mass relation for the brown dwarfs 
of (T— Orionis as determined from the observations. The purely 
theoretical 5Myr relations of Baraffe et al. (1998) and Burrows 
et al. (1997) are also shown (NEXTGEN - solid line, Burrows - 
dashed line). Note that for each calculation, the bottom point is 
an upper limit. The location of the LF dip is again highlighted 
(dotted lines). 

ing once again at / 19.5. The 50Myr NEXTGEN model of 
Barraffe et al. (1998) indicates that at the cluster distance 
of 155 pc (Barrado y Navascues et al. 2001) a member of 
magnitude / « 18 has Tefi « 2750K. This is consistent with 
the theoretical Tefi of the top edge of the sparsely popu- 
lated region in the Pleiades sequence. The /—band LF of 
the 90Myr old a-Per cluster similarly shows a dip centered 
around Mj^ ~ 12.5 (Barrado y Navascues, private comm.; 
Stauffer et al. 1999), while Figure 9 of Luhman (1999), a 
luminosity-Teff plot, hints at a sparsity of objects between 
spectral types M7-M8 in IC348. Reid & Cruz (2002) have 
derived colour-absolute magnitude diagrams for a sample of 
field dwarfs lying within 20pc of the Sun. Examination of 
their Figure 11, the M/,/ — J plot, also reveals a sparsity 
of objects in the range 2.6 < I — J < 3.0. It is clear that 
with the small number of late-M type stellar/substellar ob- 
jects in each of these studies, and in the IC2391 work the 
potential level of contamination and the location of com- 
pleteness limit of the data with respect to our claimed dip, 
individually these results are of low statistical significance. 
Nevertheless, when all the available evidence is taken to- 
gether we believe the case for a drop in the bolometric LF 
between M7-M8 is compelling. 



3 MASSES ACROSS AND BELOW THE GAP 

We have argued above for the existence of a dip in the bolo- 
metric LF of star forming regions, young clusters and the 
field near T^s ~ 2700K. These populations span a wide 
range of ages, hence the dip covers a large range of masses, 
arguing against it being related to fine structure in the IMFs. 
Instead, the obvious way to explain this is a sharp fall in 
the luminosity-mass relation between spectral types M7-M8. 
The LF can be written as in Equation 1, where dN/dm is the 



© 1994 RAS, MNRAS 000, §-| 



4 P. D. Dohbie et al. 



MF and dM/dm is the slope of the luminosity-mass relation 
where M is the bolometric magnitude and m the mass. 

dN ,dM , , 

""""^ = n^i^ 

Logically, if dM/dm increases, the LF will drop. In- 
deed, we believe that what has recently been interpreted as 
structure in the Trapezium IMF (Meunch et al. 2002) in 
fact arises from the feature in the luminosity-mass relation 
reported here. An qualitative estimate of the form of the 
7— band magnitude-mass relation for the Pleiades can be 
derived from the observations. For example, we start with 
the magnitude-mass relation from the 125Myr NEXTGEN 
model of Baraffe et al. (1998) and assume this holds true as 
far as the top of the gap in the sequence. We could equally 
use the models of Burrows et al. (1997) or D'Antona & Mazz- 
itelli (1997), but for reasons outlined in our previous work 
(e.g Jameson et al. 2002) we prefer the calculation of the 
Lyon group. Equation 2 can be readily derived by integrat- 
ing the MF from a mass mi to a lower mass m2 and again 
from m2 to an even lower mass ms, where a is the index 
of a powerlaw model MF and Ni2,N23 the number of brown 
dwarfs in the mass intervals mi to m2 and m2 to ma respec- 
tively. 

N23/N12 = {ml~°' - mj"")/(mi"° - m\~") (2) 

The ratio N23/N12 is determined from observation and 
mi and m2 are taken from the NEXTGEN model just above 
the dip. Assuming a MF for the Pleiades, dN/dm oc m""", 
a ~ 0.4, as recently derived by Jameson et al. (2002) over 
a broad mass range (O.3M0 >m> O.OSSM©), Equation 2 
can be solved for ms. Thus masses can be estimated down 
through the dip region and below as a function of I, and 
hence the magnitude-mass relation derived. We find our es- 
timates to be relatively insensitive to the assumed index a 
unless ma <Cm2. The result for the Pleiades, using only those 
candidate members with photometry consistent with them 
being single objects, is shown in Figure 3 (open circles). 
Since the component masses of the unresolved binary candi- 
dates are not known, a robust treatment of these is difficult. 
However, if we assume all of these to be equal mass binaries 
and the binary fraction above the gap, where the data are 
incomplete, to be the same as below, then the magnitude- 
mass relation is found to be nearly identical to that obtained 
by considering only the single brown dwarfs, as above. We 
choose not to proceed beyond I — 19.5 as the survey is com- 
plete over only a small area of sky at fainter magnitudes and 
the statistics are poor. The overall result, as illustrated in 
Figure 3 by the dotted line, is that beyond the beginning 
of the gap any brown dwarf mass may be greater than pre- 
dicted by either the NEXTGEN or the DUSTY models. 

This procedure has been repeated for the cr— Orionis 
cluster, using the 7— band LF data reported in Bejar et al. 
(2001). Following our argument above and since the single 
and binary sequences are not well resolved in the Bejar et al. 
work, we have ignored the influence of unresolved binarity, 
treating all their candidate members as single. As the mod- 
els of Burrows et al. (1997) are widely used to predict the 
magnitudes of planetary mass objects we opt here to use, in 
turn, both the 5Myr NEXTGEN model and the 5Myr evo- 
lutionary track of the Arizona group as the starting points 



for this calculation. The latter model has been transformed 
onto the observational plane using the bolometric correc- 
tions and Toft-colour relation of the former. As shown in 
Figure 4 we have derived the magnitude-mass relation down 
to 7 = 24.0, although due to the completeness limit of the 
survey, the point below 7 = 22 provides only an upper limit 
on mass. Except at the lowest masses where the models and 
our semi-empirical estimates appear to be reasonably con- 
sistent, we find once again that beyond the beginning of 
the gap the mass of any brown dwarf may be greater than 
predicted by the theoretical calculations. 



4 DISCUSSION 

Careful scrutiny of the evidence cited above suggests that 
the width of the LF dip is greater in the younger a— 
and 0— Orionis clusters than in either the Pleiades or the 
field populations. This is consistent with our interpretation 
in which the luminosity-mass relation undergoes a sharp 
change in form between spectral types M7-M8. For a popu- 
lation older than lOOMyrs, radius is virtually independent of 
mass right across the planetary/brown dwarf mass regime. 
However, in a population with an age of only a few Myrs, 
radius decreases rapidly with mass. Therefore, as one moves 
to lower masses and TcS s, luminosity decreases more rapidly 
in younger populations, and one should find a wider dip in 
the LF, representative of the dip spanning a constant range 
in Toff. From Figure 11 in Burrows et al. (1997) we estimate 
that the dip in the LF should be ~ 2 — 2.5x wider (in magni- 
tudes) in a 5 Myr population than in a 125Myr population, 
in agreement with the observations. 

We now speculate on the possible cause of a sharp local 
change in the form of the luminosity-mass relation in this 
Toff regime. We note that a dip in the LF observed in the 
sequences of both open clusters and the field at Mv ~ 7 (the 
Wielen gap; e.g Wielen, Jahreiss & Kruger 1983; Belikov at 
al. 1998), has been attributed to a local steepening of the 
luminosity-mass relation due to strong H~ opacity in the 
atmospheres of objects in this Toff regime (Kroupa, Tout & 
Gilmore 1990). Theoretical cool model atmospheres predict 
that dust grains of AI2O3, Fe and MgSi04 begin to condense 
in the outermost layers of a brown dwarf (or very low mass 
star) at T~ 1800K, corresponding to a Toff consistent with 
the M7-M8 LF dip (e.g AUard et al. 2001, Tsuji 2001). In- 
deed, spectroscopic observations of late-M field dwarfs lend 
support to these predictions (Jones & Tsuji 1997). However, 
the DUSTY evolutionary models of Chabrier et al. (2000) 
which include a treatment for the formation of dust, do not 
predict any discontinuity in the luminosity-mass relation 
here. In the synthetic atmospheres of AUard et al., which 
provide the outer boundary conditions for the DUSTY cal- 
culations, dust is treated as grains with sizes in the range 
0.00625 — 0.24^m. In this case, dust does not contribute 
significantly to the opacities in the near-infrared where the 
water bands are dominant (AUard et al. 2001). However, as 
illustrated by Figure 4 of AUard et al. if the grain size is in- 
creased into the regime described by Mie scattering theory, 
the opacity provided by dust is dramatically enhanced (we 
note that the opacity of dust increases with decreasing wave- 
length). While current observations argue against grains as 
large as lOO/xm, the most detailed current cloud formation 
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models predict sizes in the range 5 — 20/im (Cooper et al. 
2002). It might well be possible to reproduce the observed 
dips in the LFs with a moderate increase in the model grain 
size. 

With the recent advent of wide-field near-infrared im- 
agers (e.g FLAMINGOS) and the ever increasing availability 
of wide-field optical imagers (e.g WFI) the results of deeper, 
larger area surveys for low mass members of rich young clus- 
ters such as NGC2547, NGC2516 and IC2391 should soon 
become available. The improved statistics of these surveys 
will confirm or otherwise the reported dip in the LF be- 
tween M7-M8 and allow a more thorough examination of its 
implications for mass determinations at Tcff jS 2700K. 



5 SUMMARY 

1. We have presented evidence for a deficit of very- low mass 
stars/brown dwarfs between spectral types M7-M8 and have 
argued that this arises from a sharp local change in the shape 
of the luminosity-mass relation. 

2. Structure near the bottom end of the Trapezium IMF 
recently proposed by Muench et al. (2002) is more likely a 
manifestation of the feature in the luminosity- mass relation 
reported here. 

3. We find that brown dwarfs with spectral types later than 
M7-M8 may have masses greater than predicted by current 
theoretical models. 

4. We speculate that the onset of dust formation in the upper 
atmosphere of objects in this Toff regime is responsible for 
the change in the shape of the luminosity-mass relation. This 
might be modelled by assuming larger grain sizes than are 
currently used in the DUSTY calculations. 
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